arXiv: 1508.0406lv2 [hep-ph] 5 Nov 2015 


Vector meson spectral function and dilepton rate in the presence of strong 
entanglement effect between the chiral and the Polyakov loop dynamics 

Chowdhury Aminul Islam, Sarbani Majumder, Munshi G. Mustafa 
Theory Division, Saha Institute of Nuclear Physics, 1/AF Bidhan Nagar, Kolkata-700064, India 

(Dated: November 6, 2015) 

In this work we have re-explored our earlier study on the vector meson spectral function and its 
spectral property in the form of dilepton rate in a two-flavour Polyakov loop extended Nambu-Jona- 
Lasinio (PNJL) model in presence of a strong entanglement between the chiral and Polyakov loop 
dynamics. The entanglement considered here is generated through the four-quark scalar type in¬ 
teraction in which the coupling strength depends on the Polyakov loop and runs with temperature 
and chemical potential. The entanglement effect is also considered for the four-quark vector type 
interaction in the same manner. We observe that the entanglement effect relatively enhances the 
color degrees of freedom due to the running of the both scalar and vector couplings. This modifies 
the vector meson spectral function and thus the spectral property such as the dilepton production 
rate in low invariant mass also gets modified. 

PACS numbers: ll.15.Tk, 12.38.Mh, 25.75.-q 


2 


I. INTRODUCTION 

Quantum chromodynamics (QCD) is a theory of strong interaction that accounts for the rich phenomenology of 
hadronic and nuclear physics. However, the theory is not yet very well understood because of its nonperturbative 
nature. Ongoing ultra-relativistic heavy-ion collision experiments at RHIC BNL and LHC CERN have indicated the 
formation of quark-gluon plasma (QGP), a deconfined state of matter, transformed from hadronic matter at very 
high temperatures and/or densities as predicted by asymptotic freedom of QCD. The QCD phase diagram is not only 
essential for understanding the phenomena in the laboratory experiments involving relativistic heavy-ion collisions 
but also for the natural phenomena such as compact stars and the early universe. 

The phase diagram of hot and/or dense system of quarks and gluons predicted by the QCD has invited a lot of 
serious theoretical investigations for last few decades. The first prototype of the QCD phase diagram was conjectured 
in [1] where it looked very simple; with the passage of time more and more investigations culminated in a very 
complicated looking phase diagram with many exotic phases [2]. Nevertheless, the interest mainly revolved around 
two phase transitions - one is the chiral phase transition and the other one is the deconfinement transition. If they do 
not coincide, exotic phases such as the constituent quark phase [3, 4] or the quarkyonic phase [5, 6] may occur. So, 
an important question on the QCD thermodynamics is whether the chiral symmetry restoration and the confinement- 
to-deconfinement transition happen simultaneously or not. We note that chiral and deconfinement transitions are 
conceptually two distinct phenomena. Though lattice QCD simulation has confirmed that these two transitions occur 
at the same temperature [7] or almost at the same temperature [8]. Whether this is a mere coincidence or some 
dynamics between the two phenomena are influencing each other is not understood yet and is matter of intense 
current research exploration. 

To understand the reason behind this coincidence a conjecture has been proposed in the article [9] through a strong 
correlation or entanglement between the chiral condensate (cr) and the Polyakov loop expectation value ($) within 
the Polyakov loop extended Nambu-Jona-Lasinio (PNJL) model. Usually, in PNJL model, there is a weak correlation 
between the chiral dynamics cr and the confinement-deconfinement dynamics $ that is in-built through the covariant 
derivative between quark and gauge fields. With this kind of weak correlation the coincidence between the chiral and 
deconfinement crossover transitions [10-15] can be described but it requires some fine-tuning of parameters, inclusion 
of the scalar type eight-quark interaction for zero chemical potential /r and the vector-type four-quark interaction for 
imaginary /r. This reveals that there may be a stronger correlation between $ and a than that in the usual PNJL 
model associated through the covariant derivative between quark and gauge fields. Also, some recent analyses [16, 17] 
of the exact renormalization-group (ERG) equation [18] suggest a strong entanglement interaction between 4> and cr 
in addition to the original entanglement through the covariant derivative. Based on this the two-flavor PNJL model 
is further generalized [9] by considering the effective four-quark scalar type interaction with the coupling strength 
that depends on the Polyakov Loop (PL) field <&. The effective vertex in turn generates entanglement interaction 
between <1> and cr. Such generalization of the PNJL model is known as Entangled-PNJL (EPNJL) model [9]. This 
EPNJL model has been used to study the location of the tricritical point at real isospin chemical potential [9] and 
on the location of the critical endpoint at real quark-number chemical potential [9, 19, 20]. It has also been used to 
study [21] the effect of dynamical generation of a repulsive vector contribution to the quark pressure. The EPNJL 
model has further been generalized to the three-flavor phase diagram [22] as a function of light- and strange-quark 
masses for both zero and imaginary quark chemical potential. 

It is well known that many of the hadron properties are encoded in the correlation function and its spectral repre¬ 
sentation. The properties of the vector current correlation function and its spectral representation in the deconfined 
phase have been studied to understand the nonperturbative effect on the vector current spectral properties, e.g., the 
dilepton production rate in lattice QCD (LQCD) framework [23]. Recently, in PNJL model that takes into account 
nonperturbative effects like chiral and confinement dynamics, we have analysed [24] the effect of isoscalar-vector in¬ 
teraction on the vector meson spectral function and its various spectral properties {viz., dilepton production rate^ 
and the quark number susceptibility (QNS) associated with the conserved density fluctuation) in a hot and dense 
medium. In present article, we consider the idea of the EPNJL model in which the effective vertex generates a strong 
entanglement interaction between the chiral condensate a and the Polyakov loop <I> to re-explore the vector spectral 
function and the spectral property such as the dilepton production rate previously studied in [24]. Because of this 
strong entanglement between $ and cr, the coupling strengths run with the temperature and chemical potential. First 
we study the characteristics of mean helds with various constraints: with and without the isoscalar-vector interaction 
in both PNJL and EPNJL models. Then we further demonstrate the effect of the entanglement on vector meson 
spectral function and dilepton rate. 
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We also note that both the dilepton and real photon rate have been computed in a matrix model of QGP by considering only the 
confinement effect [25, 26] and taking into account both the confinement and chiral symmetry breaking effects [27]. 
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The paper is organized as follows: in sec. II we briefly outline the usual PNJL model and extend it with the 
entanglement effect, namely the EPNJL model. In sec. Ill we write the expression for the vector spectral function 
and its various spectral properties following our earlier calculation in Ref. [24]. In sec IV we discuss our results and 
finally we conclude in sec. V. 


II. EFFECTIVE QCD MODEL 
A. PNJL Model 

We start with the two flavour PNJL model Lagrangian with isoscalar-vector interaction [24] 


>CpNjL = '0(0 - -m-o -f 7 om) 0 + -^[('00)^ + (0i75'?0)^]-^(07^0)^ 

-^^($[A],$[A],T), 


( 1 ) 


where 0 denotes the two flavour quark field, mg =diag(mu, m^) with = rrid and r’s are Pauli matrices. = 
— igA'^\a/‘2., A'^ being the SU{3) background fields, Aa’s are the Gell-Mann matrices and g is the gauge coupling; 
Gs and Gy are, respectively, the coupling constants of local scalar type four-quark interaction and isoscalar-vector 
interaction, which do not run. U is the Polyakov potential that depends on the Polyakov Loop and it’s charge 
conjugate 3>. 

The corresponding thermodynamic potential is obtained as 


0NJL = 4>, T) -I- - 


- 2NfT 

- 2NfT 


(27r)3 

d^p 


14- 3 -f 

d^p 


- Kr4ln[J($,$)] -2V/V, 


(27r) 


3 0 ■ 


g (Ep g 3(Ep /i)/T' 

g —(Ep+/i)/T _|_ 3(Ep+/j,)/T 


( 2 ) 


Here Ep = -|- M'j is the energy of a quark with flavor / having constituent mass or the dynamical mass Mf and 

A is a finite three momentum cut-off. This Mf and the effective quark chemical potential fl are related with the scalar 
(a) and vector (n) condensates as 


Mf =mo- Gsa, (3) 

and 

0 = M - Gvn, (4) 

respectively. It is noteworthy that the value of Gs along with the current quark mass (toq) and the three momentum 
cutoff (A) are fixed in the NJL model itself to reproduce some zero temperature results namely pion mass, pion decay 
constant and the quark condensate [14]. But the value ^ of Gy cannot be fixed within the formalism of this model, 
since its value is to be fixed using the mass of the p meson which is beyond the maximum energy scale A of the model. 
So we keep Gy as free parameter and consider different choices as Gy = x x Gs, where a; is a multiplicative factor 
chosen in the range from 0 to 1. We use U from reference [14] which is fitted to lattice QCD in pure gauge theory at 
finite temperature and is given by 


G($,l»,r) b2(r). 


h. 


( 5 ) 


^ Some efforts have also been made to estimate the value of Gy mainly by fitting lattice data through two-phase model, which is not very 
conclusive. Interested readers are referred to Refs. [20, 28] and references therein. 













4 


with 

62 (T) = ao + ai(^)+a2(^) + ^3 (^ y ) ' 


( 6 ) 


Values of different coefficients uq, oi, 02 , 03 , 63 , 64 and k have been tabulated in [24]. The Vandermonde term 
./(d), $) is given as [13] 


J[$,$] 



+ 4(4)^ + 4)^) - 3(4)4')^ 


( 7 ) 


In the pure gauge theory the Polyakov potential is fitted to lattice QCD that yields a first order phase transition 
at To = 270 MeV. With this value of Tq for zero chemical potential we get, for 2-flavour case, almost a coincidence 
between the chiral and deconfinement transitions^ (Tg- = 233 MeV and T^ = 228 MeV). So the two transitions 
almost coincide [e.g.. Fig. 2(a)] but at a value higher than the range provided by the 2-flavour^ lattice QCD [29, 30] 
which is To- ~ T$ (173 ± 8 ) MeV. In Ref [14] the value of Tg was changed to 190 MeV but keeping all the other 
parameters same and obtained a lower value of (~ 200 MeV) and T$ (« 170 MeV) [e.g., Fig. 2(b)]. Taking the 
average of the two while defining gives a value almost within the range provided by the lattice QCD but then the 
coincidence is lost. In this article we work with the same Polyakov potential but the entanglement between the chiral 
and deconfinement mechanism is introduced in the next subsec. IIB. 


B. EPNJL 

The PNJL model has a weak correlation between the chiral (tr) and the deconfinement (4> and 4>) dynamics through 
the covariant derivative between the quark and the gauge fields. In addition to this there may be a strong entanglement 
interaction between d ’ and <7 as suggested by some recent analyses [16, 17] of the ERG equation [18]. Based on this 
the two-flavor PNJL model is further generalized by considering the effective scalar [9] and vector [20] type four-quark 
interaction with the coupling strengths that depend on the Polyakov field 4>. The Lagrangian in EPNJL will be the 
same as that in ( 1 ) except that now the coupling constants Gs and Gy will be replaced by the effective ones Gs)^) 
and Gv(4>). The effective vertices Gs{^) and Gv(4>) in turn generates entanglement interaction between d’ and a 
and their forms are chosen [9, 20] to preserve chiral and symmetry as given by 

Gs(4>) = Gs[l - ai4>4> - a 2 (^'^ + $^)], (8) 

and 

Gy(4>) = Gy[l - ai4>$ - q;2(^'^ + $^)]- (9) 

We note that for ai = a 2 = 0, G 5 ( 4 >) = Gs and Gv(4)) = Gy, the EPNJL model reduces to PNJL model. Also 
atT = 0, 4> = $ = 0 (confined phase), then Gs = Gs and Gy = Gy. Due to the reason already mentioned in 
the previous subsection, here again the strength of the vector interaction is taken in terms of the value of Gs as 
Gy = X X Gs, which on using (9) reduces to 

Gy(4>) = a; X Gs[l - ai4>4> - a 2 (d’^ + ^^)] = xx Gs(4>). (10) 

Now in EPNJL model, ai and are two new parameters, which are to be fixed from the lattice QCD data. The 
thermodynamic potential Depnjl in EPNJL model can be obtained from (2) by replacing Gs with Gs(4>) and Gy 
with G;/(4>). For the EPNJL model we take same values of the parameters as those in PNJL model [24] except 
the value of Tq, which is taken as 190 MeV. Then we fix the values of parameters a\ and 0 L 2 so as to reproduce 
the coincidence of chiral and deconfinement transitions within the range given by lattice QCD data at zero chemical 
potential [29, 30] and it is found that (ai, 02 )= (0.1, 0.1). We further mention that the coincidence of T„ and T <5 are 
preserved [e.g.. Fig. 2(c)] within the parameter region ai, a 2 ~ 0.10 ± 0.05. Note that the values ai and a 2 in our 
model differ from that of Ref. [9] because of the choice of different Polyakov loop potential. We chose the form of the 
potential as given in Refs. [14] whereas that used in Ref. [9] is taken from Ref. [34]. It is also noteworthy that the 
two forms of Polyakov Loop potentials are consistent with each other in the validity domain of the model [35] . 


® We note that the chiral transition temperature To- is obtained from the peak position of the da/dT whereas the deconfinement transition 
temperature T$ is that from the d^/dT. 

^ It is worth mentioning here that the chiral transition temperature is found to be Tc = (154 ± 9) MeV in the recent (2-1-1) flavour LQCD 
computations by HotQCD collaboration [31]. In (2+1) flavor QCD the chiral order parameter contains both the light quark condensate 
and the strange quark condensate. Only the former is used to define the chiral transition temperature, as the strange condensate varies 
very smoothly [32]. Now, the behaviour of the light quark condensate in (2+1) flavour and 2-flavour QCD will be similar if the light 
quark masses are similar but will be different at quantitative level as it leads to two different chiral transition temperatures simply 
because one has two different scales in the theory. The value Tc = (154 ± 9) MeV was extracted [31] entirely in reference to the chiral 
phase transition for (2+1) flavour QCD. Further, we also note that the Wuppertal-Budapest collaboration [33] has also extracted three 
somewhat different values of Tc ranging from 147 MeV to 157 MeV, depending on the chiral observables considered for the purpose. 
Since we restrict our calculation only to 2-flavour case, we stick to the corresponding Tc = (173 + 8) MeV as extracted for 2-flavour case 
in LQCD simulation [29, 30]. 
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III. VECTOR MESON SPECTRAL FUNCTION AND SPECTRAL PROPERTIES 
A. Resummed vector meson spectral function in ring summation 


The resummed vector meson spectral function in presence of isoscalar-vector interaction within ring approximation 
is in Ref. [24] 


av{u},^ 


1 

TT 


ImC'oo(w,(^ 




( 11 ) 


where Q = the four momentum of the vector meson. The imaginary part of the temporal component of the 

resummed correlator Coo is given as 


ImCoo = 


Imlloo 


l-Gv($)(l-^)Renoo 


Gv($)(l - ^)Imnoo 


( 12 ) 


and imaginary part of the spatial vector correlator is given as 


IvaCii = 


Imllii - ^Imlli 


00 


1 + - Gi^^ReHoo 




Imllii - ^Imlloo 


)^ImGoo. 

r 


(13) 


The various expressions for one-loop self-energies, IIoo and 11^^, are explicitly computed in our earlier work in [24]. 


B. Vector Spectral Function and Dilepton Rate 


The vector meson spectral function, cry, and the differential dilepton production rate are related as 


dR 

d^xd'^Q 


5a^ 1 


1 

gi^/T _ I 


CTy(w,^) , 


where the invariant mass of the lepton pair is 


and a is the fine structure constant. 


(14) 


IV. RESULTS 
A. Mean Fields 

1. Without the isoscalar-vector interaction (Gv = 0 ^ 
The gap equation for the thermodynamic potential is 


dfl(^E)PNJL „ 


(15) 


The thermodynamic potential is minimized with respect to mean fields A; with X representing cr, <f>, <i> and n. In 
this section we compare the variations of the mean fields in PNJL model with that of EPNJL one without the effect 
of isoscalar-vector interaction i.e., Gy = 0. As discussed in subsec. IIA the scalar type four-quark coupling strength 
(G 5 ) in NJL/PNJL model is fixed along with three momentum cutoff A and bare quark mass mo to reproduce known 
zero temperature chiral physics in the hadronic sector. We note that in principle it should depend on the parameters 
T and /i but it is not usually considered in NJL model [36, 37]. However, in PNJL model the PL field ($) is related 
to the temporal gluon which should make Gs to depend on $. But this dependence is also neglected in the same 
spirit [38]. So, the value of Gs remains fixed as represented by solid line in Fig. 1(a). 

Now we pay attention to the features of EPNJL model in Fig. 1(a). As soon as one introduces the $ dependence 
in the scalar coupling strength through ( 8 ) in EPNJL model, it (G 5 ) becomes dependent on both T and /x. This 
running is due to the gap equation in (15), which is solved in a self-consistent manner for different mean fields. As 
can be seen the increase in T causes Gs to decrease for a given /x and the decrease becomes faster as one increases fi. 
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This can be understood from (8) as for a given T if one increases /r, the Polyakov loop fields (<f> and $) increase and 
thus Gs decreases. Fig. 1(b) displays the temperature dependence of the scaled constituent quark mass and PL fields 
for both PNJL and EPNJL models at ^ = 0. Here we mention that for ^ = 0, $ = = |<I)| [14]. It clearly shows a 

considerable change in the chiral condensate {a = ('0'i/')) the Polyakov loop fields in EPNJL model as compared 
to those in PNJL model. For nonzero chemical potential similar behaviour of a and 4), $ is also observed. This is 
obviously due to the running of the coupling Gs which is arising due to the entanglement effect as shown in Fig 1(a). 




FIG. 1. Variation of (a) scalar type four-quark coupling strength Gs(4>) with temperature T for different values of /r and (b) 
the constituent quark mass scaled with its zero temperature value and Polyakov Loop fields with T for /r = 0 for both PNJL 
(solid lines) and EPNJL (dotted lines) model. 

Fig. 2 displays the variations of ^ and ^ with the temperature at ^ = 0 for various model conditions as discussed 
in subsecs. IIA and IIB in details. We note that To- and coincide for EPNJL model at w 184 MeV (e.g, Fig. 2(c)), 
which is almost within the range, Tc = (173 ± 8) MeV, given by the two flavour lattice QCD [29, 30]. 





(c) 


FIG. 2. Plot of dajdT and d^jdT as function of T with fi — 0 for (a) PNJL model with To = 270 MeV in Ref. [13, 14], (b) 
PNJL model with To = 190 MeV in Ref. [14] and (c) present calculation in EPNJL model with To = 190 MeV. For details it 
is referred to text in subsecs. IIA and IIB, respectively. 

We note that once ^ is introduced in the system the transition temperatures (both chiral and deconfinement) get 
reduced, which is expected. Now for a given T and /r ^ 0, 4> $ [39] generates two separate but close values 

of inflection points leading to different and T^. In that case one can take the average of and as the 
deconfinement transition temperature. For /i = 150 MeV, we found T<i> = 166 MeV and Tj, = 160 MeV and the 
average of the them (163 MeV) is very close to the value of = 167 MeV. With the increase of fi the transition 
temperatures further get reduced; for example at /r = 200 MeV, T$ = 153 MeV, = 151 MeV and = 153 MeV. 

We now discuss the differences in quark number density in EPNJL model with that of the PNJL one. In Fig. 3 
we observe that for temperature beyond 150 MeV the quark number density rises very sharply for EPNJL model 
as compared to PNJL one. This can be understood from Pig. 1(b) in which the value of Polyakov loop field rises 
very sharply beyond T = 150 MeV for EPNJL model. This indicates that the Polyakov loop field provides a strong 
correlation among the quarks at low T whereas the strength of the correlation among the quarks decreases when the 
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FIG. 3. Comparison of variations of scaled quark number density between PNJL and EPNJL model for /r = 0.15 GeV. 


value of the Polyakov loop field increases at high T and we have more and more free quarks in the system for EPNJL 
model as compared to the PNJL one. 


2. With the isoscalar-vector interaction (Gy 


Now we deal with the same set up but the isoscalar-vector interaction (Gy) is turned on through (9). In EPNJL 
model both couplings in (8) and (9) are entangled and run with T and p, by virtue of the gap equation in (15). We 
choose three different values of the strength of the isoscalar-vector interaction to demonstrate its effects within the 
EPNJL model. These values are taken in terms of Gs and the reason for which is already mentioned in the sec. 11. 



FIG. 4. Variations of (a) scaled constituent quark mass and (b) Polyakov loop fields with temperature for three different values 
of Gy at /i = 0.15 GeV in EPNJL model. 

In Fig. 4(a) the variation of the scaled constituent quark mass is shown for p, = 150 MeV. As one increases the 
strength of the vector interaction the rate of mass variation with the temperature becomes slower. Since the couplings 
run in the EPNJL model the effect of the vector interaction is more prominent than that of the PNJL model with 
fixed values of couplings [24]. In the right panel (Fig. 4(b)) the variations of the Polyakov loop fields with temperature 
at /r = 150 MeV are shown. We observe that with the increase of the value of Gy the rate of increase of Polyakov 
loop fields with temperature decreases. The differences in the constituent quark masses or the Polyakov loop fields 
for different values of Gy are however more prominent within the temperature range 165 < r(MeV) < 210. 

We have already discussed the effects of chemical potential on the transition temperatures in the previous section. 
Here in Table I we present the variations of the transition temperatures by the inclusion of the vector interaction. It 
shows that as we increase the strength of the vector interaction for the same chemical potential, the values of and 
as well as the average of T$ and Tj increase [19]. 

In Fig. 5(a) the variation of Gs with temperature at /i = 150 MeV is shown for different value of Gy. It is found 















Values of Gv and 

T-T rri_ r/n 

^ ^ <p <T 

X = 0, p = 150 MeV 

166 160 163 167 

X = 0.2, p = 150 MeV 

167 160 163.5 170 

X = 0.5, p = 150 MeV 

168 159 163.5 173 

X = 0.8, p = 150 MeV 

169 159 164 175 


TABLE I. Values of To-, T* and Tj, for different values of Gv and fi = 150 MeV. 



T[GeV] T[GeV] 

(a) (b) 


FIG. 5. Variations of (a) scalar type four-quark coupling strength and (b) scaled number density with temperature for three 
different values of Gv at /r = 0.15 GeV in EPNJL model. 


that the value of Gs increases as the strength of the vector interaction increases for a given value of temperature and 
chemical potential. This can be understood from figure 4(b) where the values of PL fields decrease with increase of 
Gy- This in turn leads to an enhancement of Gs according to (8). In Fig. 5(b) the variation of scaled quark number 
density with temperature is displayed for same fi and Gy as in Fig. 5(a). For a given temperature and chemical 
potential the number density is found to decrease with the increase of Gy. This is because the number of free quarks 
in the system is reduced since the correlation among quarks increases due to the decrease of PL fields with the increase 
of the couplings. 


B. Vector spectral function and dilepton rate 

1. Without the isoscalar-vector interaction (Gy = Q) 


Now we will be discussing the entanglement effect on the spectral function vis-a-vis the dilepton rates without the 
inclusion of the vector interaction but considering only the scalar type interaction. In Fig. 6(a) the spectral functions 
with zero external momentum [q) and zero chemical potential (^) i.e., q = fj. = 0, for PNJL and EPNJL model along 
with the free case are displayed whereas those in Fig. 6(b) are for g = /i = 200 MeV. The corresponding dilepton rates 
are shown in Fig. 7. Due to the entanglement effect through scalar type interaction the spectral function vis-a-vis 
dilepton rate for EPNJL model gets suppressed compared to PNJL model but is still higher than the Born rate. This 
could be understood in the following way. Usually the color degrees of freedom are suppressed in PNJL model due 
to the nonperturbative effect of the Polyakov Loop that causes an enhancement [24] of the dilepton rate compared 
to the Born one. As soon as the entanglement effect is introduced through the scalar type interaction that relatively 
enhances the color degrees of freedom in the system due to the running in Gs as evident from Fig. 3, hence the 
dilepton rate is reduced compared to that in PNJL model. 
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FIG. 6. Scaled spectral function in PNJL and EPNJL model are compared with the free case for (a.) fi — 0 and the three 
momentum q = 0 and (b) /j, = q = 0.2 GeV at T — 0.25 GeV with a; = 0. 



FIG. 7. Dilepton rate as a function of M/T at T = 0.25 GeV for PNJL and EPNJL model at (a) /r = g = 0 and (b) /j, = q = 0.2 
GeV with X = 0. The leading order perturbative dilepton (Born) rate is also shown. 


2. With the isoscalar-vector interaction (Gv 0) 

The free spectral function, in general, has a peak that appears at infinite value of M. This is also true for four-quark 
scalar type interaction as seen above. However, in presence of isoscalar-vector interaction Gy the peak appears at 
finite M in the resuinmed spectral function in (11) for given Gy and T: (a) below the kinematic threshold, M < 2Mf, 
the resummed spectral function has a (5-like peak due to the pole that can lead to bound state information of the vector 
meson and (b) above the threshold M > 2Mf the resummed spectral function picks up a continuous contribution 
along with a somewhat broader peak®. We here concentrate on the continuous contribution (M > 2Mf) of the spectral 
function above Tc that provides a finite width to a vector meson which decays to lepton pairs. Now we focus on the 
effects of entanglement on spectral function and dilepton rate when the vector interaction is included in addition 
to the scalar type interaction. In the left panel (Fig. 8(a)) the scaled spectral functions at ^ = 200 MeV in PNJL 
(solid line) and EPNJL (dotted line) model are shown for three different values of T. In Fig. 8(a) the peak of the 
vector spectral function, for a given T and Gy, is found to be suppressed and shifted to a higher M in EPNJL model 
compared to PNJL one. This is purely due to the entangled vector interaction as the correlation among the quarks 
in the deconfined states becomes weaker in EPNJL model. In particular, the suppression is larger at lower value of T 
and becomes smaller with the increase of T. 

The right panel (Fig. 8(b)) displays the same quantity for three different values of ^ but at a given T = 250 
MeV. Comparison with the left panel reveals that the variation of the suppression of the spectral function due to 


® The width of the peak will depend on the value of T. If T is around Tc the peak will still be sharp around the threshold [24] . 
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FIG. 8. Comparison between scaled spectral functions in PNJL (solid lines) and EPNJL (dotted lines) model for (a) a given 
chemical potential but different temperatures and (b) a given temperature but different chemical potentials at a; = 0.5. 


entanglement is strongly temperature dependent than the chemical potential. 




M/T 


(b) 


FIG. 9. Plot of (a) scaled spectral function and (b) dilepton rate as a function of M/T for PNJL (solid line) and EPNJL 
(dotted line) model for T — 0.25 GeV and p = 0.2 GeV at three differeirt choices of Gv- 


In the left panel (Fig. 9(a)) the spectral functions in PNJL (solid line) and EPNJL (dotted line) model at T = 250 
MeV and /i = 200 MeV for three different choices of Gv are compared. As evident for any value of Gy the strength 
of the spectral function for PNJL model is greater than that in the EPNJL one. In EPNJL model both couplings 
are strongly entangled through the mean fields and as one increases the strength of the vector interaction (Gy) that 
enhances the strength of the both running couplings. This in turn provides an enhancement in the strength of the 
spectral function that decays to the dilepton pairs in the medium. The entanglement effect becomes more prominent 
than that with only the scalar type interaction. These features are well reflected in the right panel (Pig. 9(b)) where 
the corresponding dilepton rates in PNJL and EPNJL models are compared. For a given Gy there is more lepton 
pairs at low mass in both EPNJL and PNJL model compared to the leading order (Born) rate. Moreover, EPNJL 
model produces less lepton pairs than PNJL one. This is due to the entangled vector interaction that reduces the 
correlation among the quarks in the medium. However, as the strength of the vector interaction increases, there is a 
relatively more dilepton production in both models. 


V. CONCLUSIONS 

In general PNJL model contains nonperturbative information of conhnement/deconfinement dynamics through the 
Polyakov Loop fields in addition to the chiral symmetry breaking dynamics. This model also employs the coupling 
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of local scalar type four-quark interaction as well isoscalar-vector interaction. The scalar type four-quark coupling 
strength is fixed along with three momentum cutoff A and bare quark mass mo to reproduce known zero temperature 
chiral physics in the hadronic sector. However, the value of the vector coupling is difficult to fix from the mass 
scale (p-meson mass) which is greater than the intrinsic scale A of the effective theory. Some efforts have been made 
in the literature to estimate the value of the vector coupling mainly by fitting lattice data. However, there exists 
ambiguity about its value as discussed. Nevertheless, the introduction of vector interaction in heavy-ion physics is 
important for study of the spectral property like dilepton rate at non-zero chemical potential. On the other hand, 
in nuclear astrophysics the formation of stars with quark matter core depends strongly on the existence of a quark 
vector repulsion. However, in PNJL model both the couplings are considered to be constant in the literature. Since, 
this model also contains temporal gluons these couplings, in principle, should depend on the Polyakov Loop helds. 
But this dependence is usually neglected and the correlation of the Polyakov Loop and chiral dynamics is a weak one 
as it arises through the covariant derivative that couples the quark and the temporal gauge field in the model. 

In this article we have extended the usual PNJL model by introducing a strong entanglement between the chiral (cr) 
and the Polyakov Loop dynamics ($), known as EPNJL model in the literature. The strong entanglement has been 
introduced via effective four-quark scalar type interaction that obeys the centre symmetry, Z{3) of pure SU{3) gauge 
group. Since the Polyakov Loop and chiral fields run with temperature and chemical potential, the entanglement 
makes also those coupling run. This entanglement effect is capable of reproducing the coincidence of chiral (To-) and 
deconfinement (T$) transition temperature within the range provided by the 2-flavour lattice data. 

The spectral function of the vector current-current correlation is related to the production of lepton pairs, which 
is considered as an important probe of the deconfined hadronic matter and has been measured in high energy heavy- 
ion experiments [40, 41]. On the other hand, at RHIC and LHC energies the maximum temperature reached of a 
hot and dense strongly interacting matter created is not very far from the phase transition temperature Tc and is 
nonperturbative in nature. In LQCD framework the dilepton production rate [23] at finite temperature but zero 
chemical potential has also been computed using a spectral function obtained from Euclidean correlation function 
through a probabilistic method that involves certain uncertainties and intricacies [24]. In our previous calculation [24] 
within PNJL model, the influence of the four-quark scalar and isoscalar-vector interaction without entanglement effect 
on the spectral function vis-a-vis the dilepton production was studied. In this article we have updated the spectral 
function and the dilepton production rate within the EPNJL model that takes into consideration the entanglement 
between the Polyakov loop and the chiral dynamics through scalar and vector interaction. In PNJL model both scalar 
and vector couplings do not run and the dominance of the Polyakov Loop fields substantially suppresses the color 
degrees of freedom around the phase transition temperature. On the other hand EPNJL model introduces a strong 
entanglement between the chiral and the Polyakov loop dynamics which relatively enhances color degrees of freedom 
in the deconfined phase compared to the PNJL model. Because of this the strength of the vector spectral function 
is suppressed and the peak is shifted to a higher energy compared to that of PNJL model but the strength is higher 
than the free one at low energy. Since the dilepton production is related to the vector spectral function, it is also 
suppressed in EPNJL model compared to the PNJL model but is more compared to the Born rate (leading order 
perturbative one) in the deconfined phase. This indicates relatively less production of lepton pairs at low energy 
with entangled vector interaction. However, as the strength of the vector interaction is increased there is a relative 
increase in the strength of the spectral function in both EPNJL and PNJL model, which also results in a relatively 
more production in lepton pairs at low invariant mass. 
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